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The discovery of the top quark in 1995 at Fermilab [l| 
concluded a long search which occupied particle physi- 
cists for about two decades. Nevertheless, a decade after 
the top quark discovery there is still much to investigate 
about this massive but fundamental particle. There are 
several ways in which new physics can make an appear- 
ance through the top quark. The large mass in particular 
suggests that the top quark may play a special role in the 
dynamics of electroweak symmetry breaking and act as 
a powerful probe in this physics sector. 

Among the technicolor scenarios, [2| predicts the ex- 
istence of a heavy meson rjT produced in gluon-gluon 
interactions (gg — > tjt — * tt). Topcolor models were in- 
vestigated even before the actual discovery of the top [3| , 
and then extended to topcolor assisted technicolor [J]. 
These scenarios account for the spontaneous breaking of 
electroweak symmetry by introducing new strong dynam- 
ics, which would explain the large top quark mass. An 
SU(3) x SU(1) residual global symmetry below 1 TeV 
predicts the existence of new particles like "topgluons" 
and a color singlet vector particle (Z') which couples pri- 
marily to the third quark generation. The existence of a 
narrow leptophobic Z' decaying to tt has been predicted 
in [U and searched for in Tevatron Run I data ](| . 

Other tt production mechanisms can be found in uni- 
versal extra dimensions theories. For example 0] inves- 
tigated a standard model in 5+1 dimensions, where a 
series of narrow resonances decay to tt final states. 

In this Letter, we present a search for narrow reso- 
nances decaying to tt by looking for an anomalous peak in 
the tt invariant mass distribution. The search algorithm 
was developed before looking at the data and optimized 
for signal from a leptophobic Z' decay. At the Tevatron, 
the Z' resonance can be produced via qq annihilation. 
The standard model (SM) predicts that the top quark 
decays more than 99% of the time into a Wb pair and 
then the W can decay either leptonically (BR = 32.4%) 
or hadronically (BR = 67.6%) @|. 

The Collider Detector at Fermilab II (CDFII) [9] is an 
azimuthally and forward-backward symmetric apparatus 
designed to study pp collisions at y/s — 1.96 TeV at the 
Tevatron. The detector has a charged particle tracking 
system [10| in a 1.4 T magnetic field directed parallel to 
the proton beam direction. This system consists of a sil- 
icon microstrip detector, which covers the radial range 
from 1.35 cm to 28 cm, and an open-cell drift cham- 
ber in the radial range from 40 cm to 137 cm. Outside 
the tracking system are the electromagnetic and hadronic 
calorimeters, with projective segmentation in rj-tp [111 ]. A 



set of drift chambers and scintillation counters detects 
muons in the central region (\r)\ < 1). The beam lumi- 
nosity is determined by measuring the inelastic pp col- 
lision rate with gas Cherenkov detectors 12| located in 
the 3.7 < \-q\ < 4.7 region. 

We search for events in which one W boson decays 
hadronically and one leptonically. The partonic final 
state for this decay channel is Ivibbqtf (£=e, /i), there- 
fore the final state events have one high momentum 
charged lepton, large missing transverse energy (^t ) U | 
due to the undetected neutrino, and four jets from the 
quarks' parton shower and hadronization processes ("lep- 
ton+jets" channel). 

The data sample for this analysis has been collected 
using high-p T triggers for electrons and muons [{J. The 
offline selection requires an electron (or muon) candidate 
with Et(pt) > 20 GeV contained in the central detec- 
tor region. A jet is defined as an energy cluster in the 
calorimeter and is reconstructed using a fixed cone algo- 
rithm with a cone of radius 0.4 in 77 — <p space. Following 
the procedure described in [3] , energies and $t are cor- 
rected to account for inhomogeneities in the detector re- 
sponse and multiple pp collisions in an event. After these 
corrections we select events containing at least 4 jets with 
\r/\ < 2 and E T > 15 GeV, and with $ T > 20 GeV. 
Jets originating from the hadronization of b quarks are 
identified (6-tagged), with about 40% efficiency, by re- 
constructing their displaced vertices The tagging 
information is not used for event selection, but only to 
reduce jet combinatorics during event reconstruction. 

Standard model processes that result in the same final 
state as the resonance are backgrounds to this search. 
The dominant components are: W boson production in 
association with jets ( W+jets), tt production, multi-jet 
events where a jet is mis-identified as a charged lepton 
and missing energy is generated by jet energy mismea- 
surement, and diboson production (WW, WZ, ZZ) with 
extra jets from initial and final state radiation (ISR, 
FSR). The expected contributions in a data set corre- 
sponding to an integrated luminosity of 680 pb -1 are pre- 
sented in Table[U Standard model tt and diboson expec- 



TABLE I: Expected number of events from standard model 
processes in a dataset corresponding to an integrated luminos- 
ity of 680 pb" 1 of data, assuming the null hypothesis. The 
observed number of events on data is 450. 



Background 


number of events 


SM tt 


199 ± 24 


diboson 


14 ± 1 


multi-jet 


22 ± 10 


W(-> tvi) + 4j 


215 ± 30 



lege, University of London, London, El 4NS, England, 9 Texas Tech 
University, Lubbock, TX 79409, r IFIC(CSIC-Universitat de Valen- 
cia), 46071 Valencia, Spain, 



tations are calculated according to their theoretical cross 
sections 15j, |16[ and their acceptances estimated from 
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Monte Carlo simulations (see later). The quoted uncer- 
tainties include the uncertainties on the theoretical cross 
sections, integrated luminosity, and acceptances. Given 
the lack of precise theoretical predictions for W+jets 
and multi-jet production cross sections, the contributions 
from these backgrounds listed in Table U are obtained 
by constraining the total numbers to add up to the 450 
events observed in data. The multi-jets to W^+jets ratio 
in this sample is fixed at 10%, which is consistent with 
an estimate obtained in a previous CDF analysis Q . The 
quoted uncertainties for these two backgrounds account 
for both the Poisson statistical uncertainty on the total 
number of observed events and the uncertainties on the 
tt and diboson predictions. 

The signature for top-quark-pair resonant production 
is a peak in the invariant mass spectrum of the tt pair 
(M t t), over a smooth background from the standard 
model contributions. In the mass range investigated here, 
all of the standard model backgrounds fall smoothly with 
increasing M t i, reflecting the parton distribution func- 
tions of the incoming protons and antiprotons. 

The resolution on the tt invariant mass is limited by 
the uncertainty on the jet energy and by the unmea- 
sured longitudinal momentum of the neutrino. The ap- 
proach adopted here is to use a matrix element technique 
[17| to reconstruct M t i for an event. We integrate the 
matrix-element for standard model tt production over 
unmeasured quantities, convoluted with detector reso- 
lution functions for the jet energies. The matrix ele- 
ment implemented for the event reconstruction is the 
standard model leading order tt production and decay, 
QQ (gg) — > ti — > W + bW~b — > qq'Mveb, and it is used to 
derive the value of M t j in the same way for events in 
simulated samples and in data. 

The a priori probability density for producing a 
tt parton-level final state, 7r pa rt({pl), relative to other 
tt final states, is the normalized differential cross section 

7Tpart(M) J|d 3 K = 



1 r l 



o Jo 



dz a dz b f k (z a )fi(z b )da k i({p}, z a P, z b P), 



where fk (/;) are the proton (antiproton) parton distribu- 
tion functions, doki the tt differential cross section, P (P) 
the proton (antiproton) momentum and {p} the set of six 
final state three-momenta. Indices k, I cover the parton 
types in the proton and antiproton respectively, and a 
sum over both indices is implied. The a priori proba- 
bility density for the parton level final state {p\ and the 
measured jet quantities {j} is given by the product 

7r p artJots(M,{j}) = 

= 7r part ({p}) • T(ji\pi) ■ T(j 2 \p 2 ) ■ T(j 3 \p 3 ) ■ T(j 4 \p 4 ), 



where T(j\p) denotes the parton to jet transfer function, 
i.e. the probability that a parton of momentum p is 
measured as a jet of momentum j. Transfer functions 
are parametrized in \rj\ and momentum of the jet, and 
are derived from Monte Carlo simulated (pythia) [l8[ 
events. Different sets of transfer functions are adopted 
for 6-quarks and lighter quarks. No transfer function is 
used for the charged lepton given its negligible momen- 
tum uncertainty. Also, the measurement is not used 
to constrain the neutrino transverse momentum. 

From 7r part ,jct s (M, {j}) we build P({p}\{j}), the prob- 
ability density for the parton momenta {p}, given the ob- 
served quantities {j}. From this distribution we derive 
probability distributions for the new variable, M t j, which 
is a function of the parton level quantities {p\, by calcu- 
lating P f (x\{j}) = JPm\{j}) ; S(x- M ti m)) d{p}. 
Having an event probability distribution for M t i we 
choose the mean as the reconstructed value for that 
event, since Monte Carlo simulations showed it to be the 
best single value M t i estimator. However, given that we 
do not know which jet matches which parton we also 
sum over all possible permutations before extracting the 
mean. If the event contains jets identified as originat- 
ing from 6-quarks we sum only over permutations with 
6-tagged jets assigned to 6-quark partons. 

To produce the M t j templates we apply the reconstruc- 
tion algorithm to each signal and background sample. 
The PYTHIA 18( generator was used to simulate both 
the Z' vector resonant pro duction (M t — 175 GeV/c 2 ) 
and tt events. ALPGEN [19( was used for the simulation 
of the W boson plus parton production, with herwig 
[2fj| used to model parton showers. All generated events 
are passed through the CDF detector simulation. The 
multi-jet background template was obtained using CDF 
II data. 

To allow direct comparison with the Tevatron Run I 
results we choose a resonance width Y z* — 0.012 Mz>; 
however, the measured signal cross section upper limits 
are insensitive to width values up to Tz> ~ 0.05 Mz>- 
The algorithm shows an intrinsic resolution on M t i of 
about 25 GeV/c 2 when applied to standard model 
tt events, in case of correct partons to jet association. 
When applied to resonance samples the mass is still cor- 
rectly reconstructed in the case of proper parton-jet as- 
signment; however, incorrect parton-jet assignment leads 
to a low tail in the reconstructed mass distribution as 
explained in Fig. [TJ The spectra plotted are restricted 
to the search region (Mti > 400 GeV/c 2 ) in which the 
standard model sources fall off exponentially in M t j . 

To derive the posterior probability for o-(pp — > Z') ■ 
BR(Z' — > tt) (we will refer to it simply as "signal 
cross section") given the observed M t i spectrum, we 



build the likelihood C(n\a, v) — 



i£{bins} 



This is the prior probability of observing n, where rii 
is the number of observed events in M t i mass bin i, 
fii = a s A s T S i + ^2jNjTji is the number of expected 
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FIG. 1: Reconstructed M t j distributions for two simulated 
standard model backgrounds. We plot the samples with the 
"hardest" (tt) and the "softest" (diboson) spectrum (normal- 
ized to unity); all the others lie in between. The inset shows 
the reconstructed M t i distributions for five signal samples. 

events in the same bin, a s is the assumed signal cross 
section, A s is the acceptance, Nj is the contribution 
from the j th background and T s j, Tji is the content 
of the i th bin of the normalized signal and background 
templates, respectively. The contributions from stan- 
dard model ti and dibosons are weighted as in Table 
HI while the contributions due to FF+jets and multi- 
jets are rescaled to account for the presence of signal 
and still match the total number of events observed 

(A / w+jcts+A m ulti-jct=Adata — Asi glla l — N t f — Adiboson ) ■ If 

one uses a flat prior distribution for the signal cross sec- 
tion and integrates over the nuisance parameters v (back- 
ground contributions and signal acceptance) then Bayes' 
theorem gives the posterior probability density P(a\n). 

The posterior probability density is used to define up- 
per limits at any given confidence level (CL), together 
with the most likely value which is regarded as the mea- 
sured cross section. This procedure is repeated for 10 
resonance masses from 450 GeV/c 2 to 900 GeV/c 2 and 
95% CL upper limits are established. 

As mentioned previously, uncertainties on template 
weighting are incorporated in the prior probability by 
means of the prior densities of the nuisance parameters. 
In this way the marginalized posterior probability density 
for the signal cross section includes the acceptance and 
cross section uncertainties. However, the relative con- 
tribution of multi-jet and PF+jet events has been kept 
fixed at 10%. To evaluate the impact of the uncertainty 
on this ratio, the multi-jet component has been either 
set to zero or increased by a factor of three, in each case 
yielding negligible change on the cross section posterior 



distribution. This is because after the event selection the 
W+jet and multi-jet M t i spectra are very similar above 
400 GeV/c 2 , as shown in Fig. [TJ We assume a 10% 
acceptance uncertainty for all the resonant ti signals. 

Other sources of systematic uncertainties can affect 
both the acceptances and the templates. These "shape 
systematics" are due to the imperfect knowledge in the 
modeling of (1) jet energy scale, (2) ISR and FSR, (3) the 
Q 2 scale for W^+jets production, and (4) parton distri- 
bution functions (PDF). To account for these uncertain- 
ties we convolute the cross section posterior probability 
with a Gaussian, whose width is estimated from a new 
set of templates and acceptances corresponding to a ±1<7 
change for each systematic. The result of this convolu- 
tion smears the cross section posterior, pushing the upper 
limits to higher cross section values. For a resonance of 
mass Mz> = 700 GeV/c 2 the uncertainties on the jet en- 
ergy scale, initial and final state radiation, W — Q 2 scale 
translate into a relative increase of the expected upper 
limits of 5%, 5%, 2% and 2% respectively, for an overall 
impact of about 10%. The uncertainty due to the choice 
of the PDFs turned out to be negligible. 

In the data sample used for this analysis, 450 events 
passed our event selection requirements, and 302 of them 
are found in the search region M t ± > 400 GeV/c 2 . The 
M t i data spectrum is shown in Fig. [21 together with the 
background expectations (above 400 GeV/c 2 ) based on 
Monte Carlo studies. In order to establish an a priori 
sensitivity of the reconstruction algorithm we generated 
1000 simulated experiments in the null hypothesis and 
extracted, for each mass point, the 95% CL expected 
upper limit, defined as the median of the upper limits 
distribution. We also calculated the central la and 2a 
frequentist coverage bands. 

In the data, the posterior probability distributions for 
the signal cross section show no evidence of resonant 
ti production. The predicted and observed upper limits 
at 95% CL are shown in Fig. [31 together with the theoret- 
ical prediction for the cross section as a function of mass 
in the leptophobic Z' model (with Y z > = 0.012 M z >) 0. 
Based on this model we exclude a leptophobic topcolor 
resonance candidate with a mass of 725 GeV/c 2 or less. 

In conclusion, we performed a search for a narrow 
heavy resonance decaying into tt in the lepton+jets chan- 
nel using 680 pb _1 of CDF Run 2 data. We set upper 
limits on the production cross section times branching 
ratio at the 95% CL. For one leptophobic topcolor pro- 
duction mechanism we exclude masses up to 725 GeV/c 2 , 
extending significantly the Run I limit of 560 GeV/c 2 [f|. 
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FIG. 2: The reconstructed M t i spectrum (data) and the stan- 
dard model prediction in the search region above 400 GeV/c 2 . 
The inset shows the same distributions in logarithmic scale. 



[2] 
[3] 
[5] 
[6] 



[7] 
[8] 
[9] 
[10] 



F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 74, 
2626 (1995); S. Abachi et al. (D0 Collaboration), ibid., 
2632 (1995). 

E. Eichten and K. D. Lane, Phys. Lett. B 327, 129 
(1994). 

C.T. Hill and S.J. Parke, Phys. Rev. D 49, 4454 (1994). 

C. T. Hill, Phys. Lett. B 345, 483 (1995). 

R.M. Harris, C.T . Hill, and S.J. Parke, 
|arXiv:hep-ph /991 1288| (1999). 

A. Affolder et al. (CDF Collaboration), Phys. Rev. Lett. 
85, 2062 (2000); V. M. Abazov et al. (D0 Collaboration), 
Phys. Rev. Lett. 92, 221801 (2004). 

G. Burdman, B. A. Dobrescu, and E. Ponton, Phys. Rev. 
D 74, 075008 (2006). 

W.-M. F. Yao et al. (Particle Data Group), J. Phys. G 
Nucl. Part. Phys. 33, 1 (2006). 

D. Acosta et al. (CDF Collaboration), Phys. Rev. D 72, 
052003 (2005). 

A. Sill (CDF Collaboration), Nucl. Instrum. Methods A 



Median 
Central 68% coverage band 
Central 95% coverage band 
Observed 95% CL upper limits: 
— • — data 



Leptophobic Z', r z .=1 .2% M z , 



■Q g _. ,j ; Simulated 95% CL upper limits: 

m 

e : 
5i 

4l 

3! 

2: 

1: 

°450 500 550 600 650 700 750 800 850 900 

M z . [GeV/c 2 ] 




FIG. 3: Predicted and experimental 95% CL upper limits us- 
ing data corresponding to 680 pb -1 of integrated luminosity, 
together with the leptophobic topcolor Z' cross section pre- 
diction. Dark and light areas define the la and 2<r ranges for 
the expected limits. The arrow marks the mass upper limit. 



447, 1 (2000). A. Affolder et al. (CDF Collaboration), 
Nucl. Instrum. Methods A 453, 84 (2000). A. Affolder 
et al. (CDF Collaboration), Nucl. Instrum. Methods A 
526, 249 (2004). 
[11] CDF uses a (z, <j>, 9) coordinate system with the z-axis 
in the direction of the proton beam; <f> and 9 are the az- 
imuthal and polar angle respectively. The pseudorapidity 
is defined as r\ = — ln(tan|), and the transverse momen- 
tum and energy as pr = p sin 9 and Et ~ E sin 9 re- 
spectively. Missing transverse energy is defined as $t = 
— E T Ai, where hi is a unit vector in the transverse 
plane that points from the beam-line to the i th calorime- 
ter tower. 

[12] D. Acosta et al, Nucl. Instrum. Methods A 494, 57 

(2002) . 

[13] A. Bhatti et al, Nucl. Instrum. Methods A 566, 375 
(2006). 

[14] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
052003 (2005). 

[15] N. Kidonakis and R. Vogt, Phys. Rev. D 68, 114014 

(2003) ; M. Cacciari et al, J. High Energy Phys. 0404, 
068 (2004). 

[16] J. M. Campbell and R. K. Ellis, Phys. Rev. D 60, 113006 
(1999). 

[17] K. Kondo, J. Phys. Soc. Jap. 57, 4126 (1988); 60, 836 
(1991); A. Abulencia et al. (CDF Collaboration), Phys. 
Rev. D 73, 092002 (2006); V. M. Abazov et al. (D0 Col- 
laboration), Phys. Rev. D 74, 092005 (2006). 

[18] T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, 
S. Mrenna, and E. Norrbin, Comput. Phys. Commun. 
135, 238 (2001). 

[19] M. L. Mangano et al, J. High Energy Phys. 0307, 001 
(2003). 

[20] G. Corcella et al , J. High Energy Phys. 0101, 010 (2001); 
|hep-ph/0210"213] (2002). 



